Abstract. Infrared to millimetre spectral energy distributions 
(SEDs) have been obtained for 41 bright ultra-luminous in- 
frared galaxies (ULIRGs). The observations were carried out 
with ISOPHOT between 10 and 200 /im and supplemented for 
16 sources with JCMT/SCUBA at 450 and 850 ^m and with 
SEST at 1.3 mm. In addition, seven sources were observed at 
1.2 and 2.2 /im with the 2.2 m telescope on Calar Alto. These 
new SEDs represent the most complete set of infrared photo- 
metric templates obtained so far on ULIRGs in the local uni- 
verse. 

I ' The SEDs peak at 60-1 00 /im and show often a quite shallow 
Rayleigh- Jeans tail. Fits with one single modified blackbody 
yield a high FIR opacity and small dust emissivity exponent 
• P <2. However, this concept leads to conflicts with several 
Q ^ other observational constraints, like the low PAH extinction 
or the extended filamentary optical morphology. A more con- 
sistent picture is obtained using several dust components with 
/? = 2, low to moderate FIR opacity and cool (50 K > T > 30 K) 
I ' to cold (30K>T>10K) temperatures. This provides evi- 
^ dence for two dust stages, the cool starburst dominated one and 
C<1 the cold cirrus-like one. The third stage with several hundred 
Kelvin warm dust is identified in the AGN dominated ULIRGs, 
^ showing up as a NIR-MIR power-law flux increase. While 
i AGNs and SBs appear indistinguishable at FIR and submm 
' wavelengths, they differ in the NIR-MIR. This suggests that 
the cool FIR emitting dust is not related to the AGN, and that 
the AGN only powers the warm and hot dust. In comparison 
O ^ with optical and MIR spectroscopy, a criterion based on the 
SED shapes and the NIR colours is established to reveal AGNs 
among ULIRGs. Also the possibility to recognise evolution- 
' ary trends among the ULIRGs via the relative amounts of cold, 
Ci ' cool and warm dust components is investigated. 
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1. Introduction 

Ultta-luminous IR galaxies (ULlRGs) emit the bulk of their 
energy in the mid- and far-infrared with Lir > 10^^ L0. Opti- 
cal and near-infrared images show a disturbed morphology and 
signatures of interaction and merging. Spectra show emission 
lines characteristic for strong starbursts, shocks (LINERs) and 
in some cases also Seyfert types (AGNs). Sanders et al. (1988a) 
suggested that the ULIRGs form an evolutionary Unk between 
starburst galaxies and quasars. The current hterature and the 
research state are excellently reviewed by Sanders & Mirabel 
(1996), and most recently by Genzel & Cesarsky (2000). 
Here we focus on the following topics: 

1) While the luminosity of an ULIRG originates from warm 
T> lOOK and cool T>30K dust, the bulk of the dust 
in normal spiral galaxies is much colder at T w 15 K. Do 
ULIRGs also have such cold cirrus-like dust, or is all the 
dust already heated up to T > 30 K? This is observationally 
related to: How does the IR SED continue beyond 100 /um 
into the submm range? For bright but less luminous IRAS 
galaxies a three component model - with cirrus, starburst 
and AGN dust - has been proposed and successfully ap- 
plied by Rowan-Robinson & Crawford (1989). 

2) What can we learn about the evolution of ULIRGs from 
such a cold dust stage? How far can the merging state be in- 
ferred from the relative amounts of cool and cold dust? For 
example, with respect to the Sanders et al. (1988a) conjec- 
ture this means: Do AGN-ULIRGs have less cold dust than 
starburst dominated SB-ULIRGs? 

3) Can the dust mass and opacity provide clues in favour of 
or against hidden AGNs? How does this compare to the 
near- and mid-IR emission of ULIRGs whose optical and 
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MIR spectra favor their power to be dominated by either 
AGNs or SBs? Observationally this relates to measuring 
the NIR-MIR SEDs, which are currentiy known only for a 
few ULIRGs. 

These topics introduce new aspects for exploring the 
ULIRGs. In this paper, based on IR to mm photometry, we 
present a sample of observational templates, and confine the 
discussion to a few basic impUcations, which may further stim- 
ulate more detailed analyses. 

2. Tlie sample 

From the Queen Mary & Westfield College IRAS galaxy cata- 
logue (QIGC, 17711 entries, Rowan-Robinson et al., 1991) we 
chose the 48 brightest sources with (1) Fgo/jm > 3 Jy and (2) 
LfIR = L40-120Atm > 10^2 L© (based on Ho = 50km/s/Mpc). 
Throughout this paper, however, we use Hp = 75 km/s/Mpc, 
thus, some of the sources reach only about 10^^^ L0. In the 
course of the ISO mission 37 out of these 48 sources could be 
observed (77 %), all objects having z < 0.2. Thus, the sample 
is representative for the nearby bright ULIRGs. In coordina- 
tion with the other ISO Guaranteed Time programmes the two 
shghtly less luminous objects NGC 6240 and Mrk 463 were 
included. Later we observed four more objects with higher red- 
shifts (0.4 > z > 0.2, Feo /im ~ 1 Jy) which were visible to the 
satelUte. 

3. Observations and data reduction 

3.1. ISOPHOT 10-200 iim data 

The observations were carried out with ISOPHOT (Lemke et 
al. 1996), the photometer on board ISO (Kessler et al. 1996). 
A description of the ISOPHOT observation templates (AOTs) 
is given in Laureijs et al. (2000). The MIR 10 to 25 jim obser- 
vations were done in triangular chopped mode (AOT P03) with 
a 52" aperture and a chopper throw of 60" which matched the 
compact size of the targets well. The FIR 60 to 200 /im obser- 
vations were performed with the array cameras CI 00 and C200 



U. Klaas et al.: Ultra-luminous IR galaxies 



3 



Table 1. Measured flux densities in Jy as a function of wavelength in /im. The objects are listed by IRAS name, except when 
another names is more common, then the IRAS name is given in the footnote. The statistical errors from signal processing 
are typically about 5 % at 100-150 /im, and reach in the other passbands 10-20% of the quoted fluxes depending mainly on 
the source/sky brightness contrast as seen by the apertures used. In case of uncertainties larger than 30 % 3 cr upper limits are 
listed. The apertures used are 30" x 30" at 1.2 and 2.2 /im, circular with 52" diameter at 10-25 /im, 46" x 46" at 60 and 90 /im, 
184" X 184" at 120-200 /im, and circular with 8", 15", 24" diameters at 450, 850, 1 300 /im, respectively. 



Name 


1.2 


2.2 


10 


12 


15 


25 


60 


90 


120 


150 


180 


200 


450 


850 


1300 


00199-7426 






0.057 


<0.120 


0.09 


0.25 






7.50 


5.34 


4.21 


3.00 








00262+4251* 


0.0054 


0.0055 


<0.087 


0.043 


0.10 


0.36 


3.36 


2.60 


2.90 


2.40 


2.40 


<2.10 








00406-3127 


















0.77 


0.63 


0.44 


0.32 








03068-5346 






<0.063 


0.040 


0.05 


0.15 


3.60 


3.34 


3.33 


2.78 


1.72 


1.35 








03158-1-4227 






<0.165 


0.085 


<0.2I 


0.25 


5.25 


3.88 


3.19 


1.75 


<1.79 


1.05 








03538-6432 


















1.37 


1.07 


0.85 


0.62 








04232-^1436 


0.0055 


0.0106 


<0.108 


0.065 


0.07 


0.39 


3.88 


4.14 


3.50 


2.30 


1.60 


1.10 








05189-2524 






0.600 


0.800 


1.00 


3.08 


13.02 


10.70 


10.00 


6.80 


4.36 


3.84 






<0.036 


06035-7102 






0.090 


0.130 


0.15 


0.53 


5.86 


5.05 


5.15 


2.90 


2.01 


1.35 








06206-6315 






<0.072 


0.040 


O.ll 


0.24 


4.84 


4.20 


4.57 


2.90 


1.87 


1.50 








12112-1-0305 






0.065 


0.100 


0.20 


0.44 


9.94 


8.95 


8.31 


5.87 


3.93 


2.93 








MK231° 






1.425 


2.400 


2.90 


8.66 


31.68 


27.34 


24.32 


14.74 


9.75 


6.88 








MK273' 






0.100 


0.250 


0.50 


2.07 


27.45 


23.78 


19.99 


13.10 


8.69 


7.40 








MK463'= 






0.425 


0.580 


0.80 


1.79 


2.33 


1.94 


1.60 


0.99 


<0.84 


0.45 








14348-1447 






<0.096 


0.108 


<0.27 


0.42 


8.21 


7.23 


6.66 


5.61 


3.76 


2.69 


0.21 


0.024 


<0.009 


14378-3651 






<0.276 


<0.399 


<0.36 


0.41 


8.53 


7.49 


5.90 


3.60 


2.80 


2.45 






<0.010 


15245-1-1019 






<0.084 


0.051 


<0.18 


0.22 


5.25 


5.43 


5.58 


3.60 


1.94 


1.31 








15250-H3609 






<0.267 


0.238 


0.56 


1.31 


8.63 


6.66 


4.16 


2.37 


1.60 


1.16 








Arp220'* 






0.147 


0.600 


1.14 


8.28 


113.35 


111.50 


109.00 


87.89 


63.95 


54.81 








15462-0450* 






0.064 


0.100 


0.15 


0.38 


2.81 






2.80 


2.44 


<2.5 








16090-0139 






<0.252 


0.100 


<0.21 


0.20 


6.70 


5.80 


4.00 


2.75 


2.00 


1.21 


<0.13 


<0.017 




NGC6240= 






0.259 


0.750 


1.00 


3.31 


23.60 


26.70 


25.90 


18.91 


12.73 


9.00 


1.00 


0.150 




17208-0014 






0.080 


0.200 


0.25 


1.32 


32.22 


31.90 


30.00 


23.00 


17.50 


12.50 


1.07 


0.155 




17463-H5806 


















0.63 


0.58 


0.37 


0.23 








18090-H0130* 


0.0200 


0.0195 


0.187 


0.300 


<0.38 


0.59 






16.63 


14.50 


10.48 


9.74 








18470-H3234 


0.0027 


0.0036 


<0.081 


0.135 


0.15 


0.44 


4.15 


3.55 


3.20 


2.45 


1.70 


1.25 








19254-7245^" 






0.123 


0.200 


0.40 


1.32 


5.57 


5.26 


4.27 


3.06 


2.30 


1.61 






0.012 


19458-H0944* 






<0.396 


<0.540 


<0.77 


0.26 




5.82 




7.61 




7.05 








20046-0623 


0.0039 


0.0040 


<0.084 


0.125 


<0.39 


0.47 


4.50 


4.71 


3.69 


3.18 


2.79 


1.40 


<0.16 


<0.026 




20087-0308 






<0.I17 


0.070 


0.10 


0.25 


5.54 


4.85 


5.96 


4.05 


2.95 


2.30 


<0.10 


<0.010 




20100-4156 






<0.090 


0.090 


0.14 


0.34 


6.53 


5.11 


4.51 


3.36 


2.09 


1.86 






<0.014 


20414-1651 






<0.168 


<0.117 


<0.18 


0.25 


6.16 


4.22 


3.60 


2.56 


1.54 


<1.35 








ES0286-199 






<0.204 


0.275 


0.50 


2.11 


13.90 


9.26 


6.73 


4.14 


2.56 


1.80 








21130-4446 






<0.090 


0.055 


<0.45 


0.15 


4.03 


3.60 


4.50 


3.94 


3.14 


2.90 






<0.023 


21504-0628 


0.0055 


0.0057 


0.040 


0.050 


0.15 


0.28 


3.33 


2.82 


2.80 


1.52 


1.26 


0.70 


<0.11 


<0.016 




22491-1808 






<0.063 


0.080 


0.11 


0.50 


6.80 


5.45 


3.55 


2.70 


1.95 


1.70 


<0.20 


0.019 




ES0148-2'' 






0.126 


0.250 


0.36 


1.55 


12.50 


11.36 


10.00 


6.40 


4.70 


3.50 








23230-6926 






0.035 


0.060 


0.10 


0.29 


4.80 


3.98 


4.10 


3.00 


2.00 


1.70 






0.010 


23365-H3604 


0.0075 


0.0065 


<0.120 


0.080 


0.15 


0.65 


7.78 


6.76 


6.75 


4.81 


3.35 


3.01 


0.17 


0.020 




23389-6139 






<0.075 


0.031 


<0.26 


0.20 


3.88 


3.30 


4.03 


2.83 


1.99 


1.52 






<0.005 


23515-2917 


















0.80 


0.64 


0.48 


0.33 









* ISO data partly distorted and uncertain; " I2540-^5708, I3428-^5608, ' 13536-1-1836, I5327-^2340, " 16504-^0228, ^ 'Super-Antennae', » 20551-4250, 23128-5919 



in the "sparse map" off-on mode (AOT P37-39) with back- 
ground positions 3-4 ' north of the source. Beam sizes used 
for photometry are given in Table |l} The four faint sources at 
higher redshift were observed in mini-map mode (AOT P22) in 
a 4 X 2 raster and restricted to the filters at 120, 150, 180, and 
200 /im. 

The data were reduced using the PHT Interactive Analysis 
tool (PIAQ) (V7.3.3e) in standard processing mode, together 
with the calibration data set V4.0 (ISOPHOT Data User Man- 
ual V4.0, Laureijs et al. 1998). This includes correction for 
non-Unearity of the electronics, deglitching (removal of data 
disturbed by cosmic particle events), and correction for signal 
dependence on the reset interval time. To handle the signal tran- 
sients in the time series of bright FCS iUuminations, only the 

' PIA is a joint development by the ESA Astrophysics Division 
and the ISOPHOT consortium led by the Max-Planck-Institut fiir As- 
tronomic, Heidelberg. 



last half was taken, when the values approached the final signal 
level. The calibration of detector responsivity and its changes 
was performed using associated measurements of the thermal 
fine calibration source (FCS) on board. 

For the chopped measurements with the PI and P2 detec- 
tors at wavelengths 10-25 /im we inspected the sequence of 
chopper plateaux and removed outliers, with the criterion that 
the uncertainty was larger than the average three sigma of the 
measurements. By this procedure typically the first chopper 
plateau of a series and those with residual glitches were re- 
moved. 

Finally the fluxes were corrected for aperture/beam size ef- 
fects. For the 3x3 pixel CI 00 array (60 and 90 /im) the fluxes 
were derived with two methods using (1) only the central pixel 
(46" X 46") and (2) the whole array (138" x 138"). Both meth- 
ods yield essentially the same flux (with some higher noise for 
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the whole array), providing evidence that the objects are point- 
like and not resolved in the FIR. 

The error propagation in ISOPHOT data reduction is de- 
scribed in Laureijs & Klaas (1999). The statistical errors de- 
rived from signal processing are about 5-20 %, depending on 
the wavelength range and object brightness (see caption to 
Tab. |l]), but systematic errors due to absolute calibration accu- 
racy are estimated to be 30 % (Klaas et al. 1998b). To account 
for the overall uncertainty in the signal derivation as well as 
relative and absolute photometric calibration we have adopted 
a general photometric uncertainty of 30 %. 

3.2. SCUBA 450 and 850 fj.m data 

The 450 and 850 /im observations were obtained on July 1st 
and 3rd, 1999, using the Submillimetre Common User Bolome- 
ter Array (SCUBA, Holland et al. 1999) at the James Clerk 
Maxwell Telescope (JCMT) on Mauna Kea, Hawaii. In gen- 
eral, photometry with 1' chopper throw, mini-jiggle and beam 
switching was performed. We used the SCUBA narrow band 
filters in order to minimize possible CO line contributions 
(see also sect. 4.1). The atmospheric transmission was deter- 
mined every hour from measurement series at various zenith 
distances (sky dips). The conditions were excellent and stable 
(r850/.m = 0.134±0.002 on July 1st and 0.209±0.003 on July 
3rd, respectively). was extrapolated from the 850 /im 

skydips (r45o^m =0.65±0.01 on July 1st and 1.09±0.02 on 
July 3rd, respectively). Mars and Uranus served as standard cal- 
ibrators, and repeated photometry measurements gave a repro- 
ducibility of better than 5 % and 14 % at 850 /im and 450 /im, 
respectively. The observing time per source was 10 - 20min 
(actually the dome carousel driver was out of order, so that 
the targets could only be observed during their oblique rising 
passage across the dome slit pointing towards south-east). The 
data were reduced using the SCUBA User Reduction Facility 
(SURF) with special emphasis on identifying noisy bolometer 
pixels and removing the sky noise determined from off-source 
pixels. 

In order to check for possible extended submm flux, jig- 
gle maps covering 2' were obtained for three sources. These 
sources were NGC6240 and 17208-0014, both reveaUng a 
submm flux excess (as discussed below), and 20046-0623 pro- 
viding the same good upper limits for different beams. None 
of the sources seems to be extended at the spatial resolution of 
about 8" and 15" FWHM, respectively. These test cases sug- 
gest that for the other sources the bulk of submm flux is also 
contained within our photometry beams and that some possible 
extended submm flux is small and lies within the total adopted 
uncertainties of about 20-30%, in particular since NGC6240 
and 17208-0014 lie in the low redshift range of our sample. 



3.3. SEST 1300 fim data 

The 1300 /im observations were obtained between July 4th and 
9th, 1999, using the MPIfR single channel facility bolometer 
(Kreysa 1990) at the Swedish ESO Submillimetre Telescope 



(SEST) on La Silla, Chile. The system provides a beamsize 
of 24" and was used in the ON-OFF mode with dual beam 
switching and a beam separation of 68" in azimuth. Pointing 
and focus were checked at least every two hours and the point- 
ing accuracy was always better than 4". The atmospheric trans- 
mission as determined by sky dips was stable at about T1300 ~ 
0. 1 . Uranus served as a standard calibrator The observing time 
per source was about 1 hour (on-source). From repeated cali- 
brator measurements we estimate the absolute accuracy to be 
about 20-30 %. 

3.4. MAGIC 1.2 and 2.2 fitn data 

J and K-band images were obtained at the Calar Alto 2.2 m tele- 
scope on August 21 and 22, 2000, using the 256 x 256 pixel 
NIR camera MAGIC (Herbst et al. 1993) with a pixel scale of 
0'.'6416. The observations were performed in standard dither- 
ing mode with a total exposure time per source of about 10 
and 30 minutes in J and K, respectively. The seeing was about 
2" preventing the identification of spatial details in the sources 
(see Appendix A). But the conditions were photometric, in par- 
ticular during the second night. The UKIRT standard stars FS2, 
FS27 and FS35 were observed for flux calibration. The data re- 
duction followed the standard procedures. Based on the cross 
calibration of the standard stars, we estimate the photometric 
accuracy to be about 15%. In the case of IRAS 18090H-0130 
the uncertainty is somewhat higher due to superposition of two 
stars. 

4. Results 

4.1. Spectral energy distributions 

The fluxes from our observations are listed in Table[ll The 60 
and 90 fim values agree within 20 % with the IRAS 60 and 100 
fim ones. For many sources, where IRAS provided only upper 
limits at 12 or 25 /im, now the 10-25 /im fluxes could be mea- 
sured. Also, for 16 sources out of 41, submm/mm fluxes and 
good upper limits are provided, as well as NIR 1.2 and 2. 2 /im 
fluxes for seven sources. 

The spectral energy distributions (SEDs, as measured, not 
corrected for redshift) are shown in Fig. U supplemented by 
literature data. The remarkable features of the SEDs are: 

1) For each galaxy the maximum of the SED can now be 
clearly determined. It lies between 60 and 100 /tm. The 
ISOPHOT long wavelength filters beyond 100/im clearly 
outline the beginning of the Rayleigh-Jeans branch. In 
some cases (e.g. NGC6240, 17208-0014) a somewhat 
plateau-like broad maximum is revealed, suggesting a high 
opacity even in the FIR and/or the presence of various cool 
to cold dust components (see also Fig.||). 

2) For 20 sources the shape of the Rayleigh-Jeans branch 
can now be determined by the mm and submm data 

within the 



points. As discussed in Sections 3.1 and 3.2 



error budget the effects of different instrumental beams 
and a possible contribution of extended flux appear to be 
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Fig.l. 
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LINER 



z = 0.0245 
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LINER 
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Fig. 1. continued. 
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Fig. 1. continued. Spectral energy distributions of ultra-luminous IR galaxies, ordered along the columns by RA. Upper limits 
are indicated by a downward arrow and in the case of ISO, SCUBA and SEST measurements in addition by open symbols. 
The redshifts and, if known, the optical spectral types are listed. The wavelength and frequency ranges are as observed and not 
corrected with regard to the rest frame of the objects. The dotted hues represent modified blackbody fits. Emissivity exponent 
/3 and opacity rioo/im are listed in the plots. The resulting temperatures T are corrected for redshift. For those sources with no 
submm/mm data available we used an average (3= {1.6). The solid lines indicate the shapes of the SEDs, dash-dotted lines are 
used where they are uncertain. At NIR-MIR wavelengths the SED shape is obtained by connecting the data points, and at long 
FIR and submm wavelengths it is determined by the modified blackbody fit. Luminosities were derived by integrating the SEDs 
between 10 and 1000 /im. Dashed lines in the mm range indicate synchrotron extrapolations from cm observations. The literature 
data are from Sanders et al. (1988a,b), Crawford et al. (1996), Marx et al. (1994), Murphy et al. (1996), Rigopoulou et al. (1996), 
Due et al. (1997), Rigopoulou et al. (1999) and Lisenfeld et al. (2000). For some sources the NIR fluxes of the central regions are 
plotted with a "x" (Diameter = I", from ScoviUe et al. (2000), except 5'.'5 for I8090-H0I30 from Hill et al. (1988)). 



negligible. In particular for those cases (e.g. NGC6240, 
I7208-00I4), where ISOPHOT akeady indicated a broad 
maximum, the submm measurements confirm the high 
flux prediction from the 100-200 /xm measurements. Also, 
where ISOPHOT indicated a steep Rayleigh- Jeans branch 
(e.g. 16090-0139, 23365H-3604, 23389-6139) the mm and 
submm fluxes and upper limits are low, demonstrating the 
consistency between the instruments. 
3) Shortward of the maximum at around 60-100 fim the SEDs 
exhibit two basic shapes: 



i) a flat NIR plateau followed by a jump-like flux increase 
at about 10 /xm (e.g. Arp220 and 121 12-H0305). In some 
cases also indications of PAH emission around 7.7 /im 
and silicate 9.7 /xm absorption features (e.g. Arp220 
and NGC6240) are recognized. 

ii) a power-law-like flux increase from the NIR to the 
MIR (e.g. Mrk463) or FIR (e.g. Mrk231). Spectral 
PAH emission or silicate absorption features may be 
present, e.g. for Mrk463, compare also with Fig. 3 in 
Rigopoulou et al. (1999), but due to the high continuum 
level they are diluted in the broad band photometry. 
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4.2. CO line and synchrotron emission contamination of 
submm/mm fluxes 

The determination of the Rayleigh Jeans branch allows a de- 
tailed analysis with respect to the dust emissivity A^*^ and 
the opacity rioo/^m, as carried out in the next section, under 
the condition that the emission is of thermal nature. There- 
fore, beforehand one has to check possible contamination of the 
submm and mm fluxes by CO lines and/or synchrotron emis- 
sion: 

1) For all our sources observed at 1300 /im with SEST the 
CO (2- 1 ) line at 230 GHz ( 1 304 ^m) moves out of the filter 
band pass due to their redshifts of z <; 0.05. 

2) The CO (3-2) line at 345 GHz (869 ^m) may contribute to 
the SCUBA 850 /^m fluxes for those sources with z < 0.04. 
Direct CO (3-2) observations are only available for Arp220 
(Mauersberger et al. 2000). Therefore, we estimated CO (3- 
2) line strengths from the CO (1-0) fluxes (Solomon 1997, 
Downes and Solomon 1998, Gao & Solomon 1999), adopt- 
ing an intensity ratio R = Ico (3-2) /Ico (i-o) = 0.9 as 
found for Arp220. The strength of the redshifted CO(3- 
2) line, however, is reduced according to the transmis- 
sion of the SCUBA 850 ^m filter which lies between 55 % 
and 30 % for our sources with z ^ 0.04. It turns out that 
only four sources have a CO (3-2) line contribution to the 
850 /im flux which exceeds 5 %: Arp220 and Mrk 273 (both 
1 1 %), Mrk 231 (22 %), and NGC 6240 (35 %). These con- 
tributions are smaller or of the order of our adopted pho- 
tometric uncertainties. Therefore we decided not to correct 
for the CO (3-2) line contribution. 

3) In 20 cases with available radio fluxes (Condon et al. 1990 
and Crawford et al. 1996) extrapolations with spectral in- 
dices between 0.5 and 1.0 from the cm range towards 
shorter wavelengths are far below (< 10%) our measured 
submm and mm fluxes (many of the extrapolations lie even 
below the range plotted in Fig.[l|). The only exception is 
Mrk463 (with two Seyfert nuclei) where the 1.3 mm flux is 
variable, thus dominated by synchrotron emission (Chini et 
al. 1989a and Marx et al. 1994.) 

Our conclusion is therefore, that, like the FIR emission, the 
bulk of the submm/mm flux is of thermal nature being emitted 
by dust. 

4.3. Dust parameters 

In order to characterise the dust emission, the SEDs are fitted 
with modified blackbodies. Such fits, however, are not unique. 
They rely largely on the mass absorption coefficient k and its 
wavelength dependence (3, both still being a matter of debate. 
Values of (3 between 1 and 2 are commonly used (e.g. Hilde- 
brand 1983). In case of a flat Rayleigh-Jeans tail the SEDs can 
also be modelled by several dust components. Since the inter- 
pretation of the dust emission as well as the derivation of the 
dust mass depend on the blackbody models used, we investi- 
gate the two main cases. They represent simplified formalisms. 
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Fig. 2. Distribution of Tioo/^m versus j3 for those sources 
with measured submm/mm fluxes. Different symbols corre- 
spond to optical spectral types as in Table 4: + Seyfert 1, x 
Seyfert 2, i LINER, ■ HII/SB, and • for not classified; the 
quasar PG0050H-124 denoted by * is also included for com- 
parison (from Haas et al. 2000a). Arrows indicate lower limit 
cases. 

each relying on implicit assumptions, and a realistic description 
probably lies between these two extremes. In the following two 
subsections the FIR-submm range is investigated, and the MIR 
part is addressed in the third subsection. 

4.3.1. Single modified blackbody 
We used the following model: 

Sx=Bx{T)-{l-e-^^),Wi\h (1) 
T\ = Tioopm • {lQQ^i.■m/\Y (2) 

The SEDs are fitted between 60 and 1300 /im with a single 
modified blackbody leaving T, (3 and Tioo^^m free (minimising 

in a grid search). Although the emissivity exponent (3 is 
a free parameter in eqn. 2, it is kept constant over the whole 
wavelength range. In order to keep the parameter space under 
control, we have decided not to introduce a A dependence of (3. 

The relation of the parameters (3, rioo^im and T with the 
SED shapes is: 

• j3 corresponds to the slope of the Rayleigh-Jeans tail, 
whether it is flat (small (3) or steep (large (3); 

• Tioo /jm determines, whether the peak plateau of the SED is 
narrow (small rioo^im) or broad (large rioo ,„«); 

• Like for normal blackbodies (Wien's displacement law), T 
is related to the wavelength, at which the SED maximum is 
located. 
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Fits were performed only for the "mm-subsample", i.e. for 
those 22 sources with submm/mm fluxes available. We also 
used the 8 cases of upper limits^ treating them formally as de- 
tections, and after the fitting procedure taking into account that 
they provide lower limits for /3. In some cases, where the ISO 
60 or 90 /im fluxes show a large error or deviate strongly from 
the IRAS 60 and lOO^m fluxes, we also included the IRAS 
data (e.g. for 16090-0139). 

The fitted parameters /?, rioo and T are listed in Table ^ 
together with x^. A visual impression of the quality of the fits is 
given in Fig. [l|. In summary, the results for the mm-subsample 
are: 

• 1.2</5<2.0 (Fig.|) 

• 0.5<rioop,„<5 (Fig.|) 

• 50K<T<70K. 

with an uncertainty of about 0.2 for (3, 25% for r and 5- 10 K for 
T. The quoted uncertainties have been estimated considering 
"bent" SEDs with the FIR fluxes reduced and the submm fluxes 
increased by their measurement uncertainty and vice versa. It 
should also be noted that for sources with significant MIR flux 
(an extreme case is Mrk463) additional uncertainties in the de- 
termination of T and/or tiqo are introduced when using the 
full 60 /im flux, which is partially due to a dust component not 
considered here. 

The quoted parameter values should be considered with 
some tolerance and their interdependence borne in mind: 

• For two SEDs with maxima located at similar wavelengths, 
T decreases with (3. Due to the steeper fall-off at long wave- 
lengths the width of the peak plateau can only be main- 
tained, if the maximum is shifted to longer wavelengths. A 
complementary pair are 19254-7245 ((3=1.2) and 20046- 
0623 (/3=1.8). 

• For two SEDs with similar (3, but largely different tiqo ^im, 
T increases with Tioo/^m- Due to the broader width, but 
similar shape on the long wavelength side, the maximum 
is shifted towards shorter wavelengths. A complementary 
pair is Mrk231 and 14348-1447. 

On the other hand, there seems to be no correlation between 
'''loo Mm and (3 as can be seen from Fig. ^. 

The dust parameters were determined via equation (1) only 
for the mm-subsample. For the remaining sources with wave- 
length coverage limited to 200 /im (IR-subsample) [3 could 
not be fitted reliably (as we found from tests with the mm- 
subsample using only the 60-200 /im fluxes). For the IR- 
subsample we kept (3 fixed using the average value /3 = 1 .6 de- 
rived from the mm-subsample. Then rioo/jm and T could be 
determined reasonably well from the 60-200 /im fluxes alone. 
(Exceptions are 00262H-4251, 15462-0450, 18090H-0130 and 
19458H-0944 which have less complete spectral coverage due 
to bad quality measurements as flagged in Table [l|. In these 
cases Tioo/iim was fixed to 6.0). The resulting values lie in the 

^ these are: 14378-3651, 16090-0139, 20046-0623, 20087-0308, 
20100-4156, 21130-4446, 21504-0628, 23389-6139 



same range as for the mm-subsample (Table ^and Fig. |l]). As a 
check, we fitted also rioo and T of the mm-subsample with 
a fixed /3 = 1 .6 using only the 60-200 /im fluxes. The results are 
basically consistent with those obtained from the longer wave- 
length coverage, except for the sources with extremely low or 
high true (3. Hence, in the discussion below we can mostly use 
the full sample, and only where (3 plays a role, we confine it to 
the mm-subsample. 

4.3.2. Multiple modified blackbodies 

As derived in the previous section, for the majority of the "mm- 
subsample" sources (11 out of 14, not having lower limits for 
/3, one of them having rioo /^m = 5, Fig. ^) it is not possible to fit 
the FIR-submm SEDs properly with one single modified black- 
body with an emissivity law of A^^, rather the superposition of 
two or more modified blackbodies is required. 

In the low opacity case Equation (1) can be approximated 

by 

Sx=B^{T)-T^<^Bx{T)-\-P (3) 
Based on Equation (3), implicitly assuming rioopm ^ 1, 
we fitted several modified blackbodies with free T and fixed (3 
= 2 to all observed SEDs, as illustrated for some examples in 
Fig. ^. They have temperatures in the range between 30 and 
50 K (cool), and between 10 and 30 K (cold). However, the 
decomposition into the various components is not unique, be- 
cause of the sparse data coverage in the sub-mm range. 

In the case of multiple blackbodies no direct conclusion 
about the opacity rioo^m can be drawn (nevertheless, in Sect. 
5.1.2 below, Tioo/jm will be constrained using CO data). The 
most realistic case might be that of several blackbodies with 
/3 w 2, and a range of opacities from low to partly high. 

4.3.3. Description of the NIR-MIR emitting dust 

The two basic SED shapes in the NIR-MIR outhned in Sect. |1 
can be formally fitted by a superposition of several warm dust 
components. For the cases with flat NIR plateau the maxi- 
mum temperatures are about 100-150K (e.g. Klaas et al. 1997, 
1998a). The power-law-like SEDs can be approximated by a 
suite of blackbodies up to the dust grain evaporation tempera- 
tures of about 1000-1 500 K (the hotter the blackbody the less 
dust mass is involved). Modelling of the continuum is ham- 
pered by the presence of strong spectral features like Y'AH 
emission and silicate absorption. Using higher spectral resolu- 
tion, Laurent et al. (2000) investigated this spectral part quan- 
titatively. 

4.4. Luminosities 

Table ^ lists the luminosities derived within various bandpasses 
in the rest frame of the objects by integrating the spectral en- 
ergy distribution as outlined by the thick solid and dash-dotted 
lines shown in Fig. |l] for the indicated wavelength ranges. On 
the Rayleigh-Jeans tail and around the SED maximum this 
comprises the single blackbody curve obtained with the eqn. 
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Fig. 3. Examples for fits with several optically thin (Tioo^im •C 1) modified blackbodies (BBs) with an emissivity law of A~^. 
Top row: superposition of two BBs; middle row: possible superposition of three BBs; bottom row: fits can also be done for those 
SEDs with no submm data: approximation by two BBs. 



(1) fit, and shortward thereof the lines connecting the data val- 
ues by linear interpolation. 

The total IR-submm luminosity Lio_iooo/im is dom- 
inated by the FIR in the wavelength range 40-150 /um, 
while the 1 50-1 000 ^m submm range plays a minor role 
(Li5o_iooo/im <0.1 •L4o_i50Mm) as Well as the 10-40/;im 



MIR range (except for Mrk463 and the z > 0.3 sources which 
are IVIIR dominant). The luminosities Ls-iooo/im extrapolated 
from the four IRAS bands (formula cf. Table 1 in Sanders & 
Mirabel 1996) typically slightly overestimate our IR-submm 
luminosity values by about 15%; nevertheless this is still a 
good agreement. 
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The MIR/FIR luminosity ratio has a median value of about 
0.3. Thus, the sample of bright nearby ULIRGs preferentially 
comprises objects with cool MlR/FlR colours (compared with 
quasars having Lmir/Lfir > 1, cf. Haas et al. 2000a). Though 
the luminosity range of the ULIRG sample spans about one 
decade, there is no trend of luminosity with optical spectral 
type or MIR/FIR colours. 

For sources without submm/mm observations available, the 
submm luminosity Lsubmm = Li5o-iooo£im is extrapolated us- 
ing the average value /3 = 1 .6 (Sect. [4.3. 1 ), while the actual 
Lsubmm depends on the actual value of /3. A check on the 
mm-subsample shows that, for the case of minimum /3 « 1.2 
or maximum /3 « 2, Lsubmm obtained using the average (3=1.6 
can deviate from the true value by factors of 2 and 0.5, respec- 
tively. 

4.5. Dust masses 

In order to derive the dust masses, we used the standard ap- 
proach based on Hildebrand (1983) and further developed by 
various authors (e.g. Chini et al. 1986, Kriigel et al. 1990): 
An-a^-p D^-S. 



M 



dust 



''850/im 



3 • a2 ■ Qe™(850/xm, P) Bs5o^miT) 
with 

• average grain size a = 0. 1 pm 

• grain density p = 3 g cm^"^ 

• distance D 

• flux Sgso/^m (in the restframe of the object) 



(4) 



• Planck function B 



(T) at restframe 850 /im, and 



• dust grain emission efficiency Qe„j(850/im, (3) where 



(5) 



Qem{A,/3) = Qo-a-{250/A}'3 
wifli Qo = 40 cm-i for A = 250 ^m. 

The first fraction in eqn. 4 is conveniently summarized yielding 

A f - ^ -D^ ■ Ss50^l.m , 
jWdust — ■ -n 77^ C") 



with 



with 

13=' 



K 



0=2 
S50fim 



{250nm/850pm} 



f}-2 



(7) 



= 0.865 cTO^g- 



which corresponds to 



Our "choice" of k is consistent with that favoured by 
Kriigel et al. (1990, see their eqn. 10), and Lisenfeld et al. 
(2000). Note that we account for (3 in the wavelength depen- 
dence of Qom and n, respectively. In the case of /3 < 2 this leads 
to a larger value of k, hence to dust masses which are smaller 
than for [3 = 2. 

Table ^ lists the dust masses derived from flux values at 
850 /xm, where the emission is certainly transparent (and no 
mass is "hidden"). For the further discussion, three different 
kinds of dust mass are determined for each galaxy: 



(1) For the assumption that the FIR emission can be modelled 
by one single blackbody with free r, j3 and T, as described 
by eqn. 1 and shown in Fig.|l|. 

For the case of several components with fixed (3 = 2 (eq n. 3) 
the temperature decomposition is not unique (Sect. 4.3.2 ), and 
the uncertainty in the derived total dust mass becomes large, 
in particular for those sources without any submm/mm data. 
Therefore Table |^ lists two estimates: 

Q = 2 

(2) The total dust mass MJ^j^^^^^j^ associated with the multiple 
"optically thin" blackbodies. 

(3) The dust mass associated with the "bulk FIR emission com- 
ponent" Ul=f^y 



While M 



d(total) 



may be considered as a maximum dust 

= 2 

mass, MJjj^pjj^^ is derived under the assumption of one single 
opticafly thin A"^ modified blackbody of T«35-40K fitted 
to the FIR 60-200 /im range. It contributes the bulk of the lu- 
minosity and may be considered as a firm lower limit for the 
dust mass, in particular if other colder dust components exist. 

8 = 2 

'^d(totai) contains contributions from cold dust components, it 

Q = 2 

is typically much larger than MJj^pjj^^ . 

For all cases, the dust mass does not show any correla- 
tion with the total, mid- or far-infrared luminosity. But the 
dust mass is quite well correlated with the submm luminosity 

Ll50-1000/jm- 



4.6. Size of dust emitting regions 

The smallest possible extent of the FIR emitting region is listed 
in Tab.^. For the case of an opaque blackbody (eqn. 1), the 
brightness radius rb is determined via 

" U • BiooM™(r) • (1 - e--ioo."0 J 
For the case of a transparent blackbody with fixed (3 = 2 
(eqn. 3), the low opacity condition tiqo ^j.m < 1 requires that 
the extent of the dust emission cannot be smaller than a 
minimum size. Such a "smallest transparent radius" can be 
roughly estimated as follows: the dust mass (eqn. 6) - as de- 
termined for the bulk FIR emitting component at T 30- 
40 K - is distributed homogeneously in a "minimum face-on 
disk" of radius r,-. Thereby the following standard conver- 
sions are used: 0.1 M© ■pc"^ ry =0.4 (Whittet 1992), and 
Tioo urn =0.006 • TV (Mathis et al.l983). Thus 



Mdust[MQ] 



1/2 



[pc] 



(9) 



41.7 -TT 

Typically r^., for rioo^m — 1, is a factor 2-5 larger than rf, 
(see Tab.^, and we will discuss further below whether a dust 
emitting region as large as that given by ir is still consistent 
with other data. 



5. Discussion 



U. Klaas et al.: Ultra-luminous IR galaxies 



13 



Table 2. Infrared luminosities, dust temperatures, dust masses and IR source sizes. 

The luminosity distance is determined as Dl = c / Ho * (z + / 2), i.e. qo = 0, with Hq = 75 km/s/Mpc. The infrared luminosities 
Lmir (3-40 /im), Lfir (40-150 fim) and Lsubmm (150-1000 /im) are determined from the SED curves shown in Fig.|l], by inte- 
grating over the indicated wavelength range shifted to the rest frame of the objects. The dust mass is derived for three models of 
the IR emission: 

1) Mj^^'^'"^ for the single modified blackbody with free T, t and /3 as described by eqn. 1 (minimising x^) and shown in Fig.[l| (/? 
values in brackets < > means that a fixed average of /3= 1.6 is used, lower limits are indicated by a > sign). 

2) M^j,pjj^ for a single modified blackbody with fixed /3 = 2 fitted by eqn. 3 to the FIR 60-200 ^im range (for details see 
sect. O). 

3) Mj^^j^j^j^ for multiple modified blackbodies with fixed /3 = 2. 

The sizes of the dust emitting regions are given for the two cases of a single blackbody as (1) brightness radius rb, and (2) as 
"smallest transparent radius" r^., at which rioo/^m = 1 (for details see sect. i.6). For comparison M(H2) is listed (using the galactic 
conversion factor M(H2) = 4.6 • Leo)- 
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* measurements partly distorted and uncertain, therefore rioo^tm has been fixed to 6.0 

+ good determination of total dust mass, data used for estimate of gas-to-dust ratio (see Sect, 
" from Solomon et al. 1997, Table2 

'' from Mirabel et al. 1990, recalculated by Gao & Solomon 1999 
^ from Sanders et al. 1991, recalculated by Gao & Solomon 1999 



5.7. Evidence for cold cirrus-like dust showed that the observed SEDs can be modelled by either one 

single tepid (50-75 K) blackbody with /3 < 2, partly inferring a 
The observations show that, though the SEDs usually peak high opacity even in the FIR, or by several cool to cold ti'anspar- 
at 60-100/im, they can stay on a plateau bef ore the on set of 
the Rayleigh- Jeans tail. The analysis of Sect. |4.3.l| and 4.3.2 
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ent dust bodies with /3 = 2. Which of these scenarios (or which 
combination) might more realistically apply for the ULIRGs? 

In the following discussion, we will check for observational 
evidence for the existence of cold cirrus-like dust in ULIRGs 
being in line with the multiple blackbody model. One prereq- 
uisiste to allow a fit of the SEDs with several blackbodies is 
that /3 « 2. Even in the case of several blackbodies the opacity 
cannot be low at all, since the derived dust mass and the size of 
the emission area exceed other constraints. In order to reduce 
the amount of dust mass and size, the opacity of some of the 
blackbodies must be increased (whereby the dust temperature 
also rises, resulting in a lower dust mass and a smaller region). 
Thus, we end up with the picture of several blackbodies with 
/3 = 2, which are partly optically thick. 

5.1.1. Evidence against the single blackbody model 

The following arguments can be brought forward against the 
concept of one single blackbody to describe the FIR-submm 
SEDs of most of our ULIRGs (i.e all of our mm-subsample 
except Arp220): 

1) In general, the optical morphology is disturbed and sug- 
gests a variety of patchy dust complexes with less dense 
regions in between. Since in such a scenario various dust 
temperatures are also expected, multiple blackbodies ap- 
pear to provide a better physical description than one single 
blackbody. 

2) For a variety of objects several authors found observational 
evidence that « 2 in the FIR-submm range (A > 100 /im): 
l.S^P ^ 2.0 was found for interstellar dust in our Galaxy 
(e.g. Mathis et al. 1983, Lagache et al. 1998), /3w2.0 for 
normal spiral galaxies (e.g. Bianchi et al. 1999) and the 
archetypal starburst galaxy M 82 (e.g. Thuma et al. 2000), 
and (3 = 2.0±0.2 for active (Markarian) galaxies (Chini et 
al. 1989b). In contrary to that, Dunne et al. (2000) found an 
average value /3 « 1.3 from their SCUBA Local Universe 
Galaxy Survey. However, one weakness of this study ap- 
pears to be the restriction to three wavelength points at 60, 
100, and 850 /im without probing the maximum and the 
start of the Rayleigh Jeans branch in sufficient detail. One 
common object between their sample and ours is Arp 220. 
While their single temperature fit provides /3 = 1 .2 our sin- 
gle modified blackbody fit to a much better sampled SED 
gives (3=1.1. Also, the spread of (3 between 1.2 and 2.0 
in our ULIRG sample would suggest a large variety of dust 
properties, like grain size distribution, among the same type 
of object with probably similar evolution histories. It ap- 
pears much more plausible to explain the SEDs by a com- 
position of dust components all with consistent emissivity 
laws of /3 « 2. 

3) We developed the following quantitative method to check 
the consistency of the derived tiqo f^m values from the sin- 
gle blackbody fit with other extinction sensitive quantities. 
One is strength of the PAH 7.7 /im feature which is ubiqui- 
tous in ULIRGs (Rigopoulou et al. 1999). There must be a 



sufficient number of UV photons to excite the R\Hs, how- 
ever, an intense starburst is not a prerequisite, as can be seen 
from the wide distribution of PAH emission in our Galaxy 
(Mattila et al. 1996) and in NGC 891 (Mattila, Lehtinen & 
Lemke 1999). Mattila, Lehtinen & Lemke (1999) found 
that the PAH distribution is similar to the one of large 
dust grains and neutral molecular clouds. It is reasonable 
to assume that the PAH carriers are mixed with other con- 
stituents of the ISM, at least on the spatial scales of kpcs 
we are looking at with the resolution of our observations. If 
the PAHs are not cospatial with the dust component emit- 
ting the bulk of the FIR emission, then there is a second 
noticable dust component per se. 

In the case of all the FIR emission being irradiated by one 
single blackbody, we would expect a mixture of the PAH 
carriers with the FIR emitting dust and a strong decrease 
of the 7.7 /im PAH strength (published by Rigopoulou et 
al. 1999, see examples in Fig.^ and Tab.||) relative to 
the submm and FIR fluxes with increasing rioo/^m- This 
decrease should follow the direction of the Ay-vector in 
Fig.||, (using the galactic extinction curve by Mathis et al. 
(1983)). However, Fig. I shows that the normalised 7.7 /im 
PAH fluxes are largely independent of rioo/^m, except for 
Arp220 which is further discussed in Haas et al. (2001). 
Further clues on the relation between PAH and FIR-submm 
emitting dust comes from the comparison of the ULIRGs 
with normal galaxies: We examined the PAH / 100 /im and 
PAH/ 850 /im flux ratios of the ULIRGs with that for nor- 
mal galaxies, for which both 850 /im and PAH data are 
available. Among the SCUBA sample of 104 galaxies ob- 
tained by Dunne et al. (2000), 20 sources were also mea- 
sured with ISOPHOT-S (by several observers) covering a 
similar area (24") as the SCUBA observations. We derived 
the PAH 7.7 fim fluxes from the ISO data archive products 
processed with OLP Version 9.0 (Laureijs et al. 2000). The 
values are listed in Tab.|3]and an example for the MIR spec- 
tra is shown in Fig.^ Note, that now we consider only the 
vertical distribution in Fig. |] and ignore the horizontal one: 
The range for the normal galaxies is indicated by the thick 
vertical bars on the left-hand side. Strikingly the normal 
galaxies populate a similar PAH/ 850 /im range (3.9 ± 2.3) 
as the "typical" ULIRGs (3.7 ± 1 .5, excluding Arp220), but 
the PAH / 100 /im range (0.03 ±0.014) Ues higher than for 
the ULIRGs (0.01 ± 0.004). This difference cannot be ex- 
plained by PAH destruction in the ULIRGs, since it should 
be reflected in the PAH / 850 /im ratio as well. Neither can 
this difference be due to extinction, because the shift of the 
normal galaxies with respect to the ULIRGs must then be 
larger at 850 /im than at 100 /im (see the length and orienta- 
tion of the extinction arrows in Fig. || The difference of the 
PAH / 100 /im ratio between normal galaxies and ULIRGs 
can be best explained by the higher dust temperatures of the 
ULIRGs of about T = 30-60 K compared with T= 15-30 K 
in normal galaxies, which increases the 100/im flux with 
respect to that at 850 /im. This behaviour also suggests that 
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Fig. 4. Examples of ISOPHOT-S 5.8-11.6^m spectra for the 
ULIRG Arp220 and the comparison galaxy NGC 23. The dot- 
ted line indicates the continuum subtracted for the estimate of 
the PAH 7.7 /im peak flux strength. Arp220 has a strong sil- 
icate 9.7 /im absorption feature (e.g. Dudley 1999), which is 
accounted for in the continuum estimate. 

the PAH carriers not only trace the dust in active regions, 
but are even more related to the entire dust content. 
4) We compared the dust mass with the 

molecular gas mass derived from integrated CO luminosi- 
ties (M(H2) in Tab.^). The resulting gas-to-dust mass ratio 
lies in the range between 500 and 2500, with an average 
of about 1500 ± 500. This is much higher than the canon- 
ical value of about 150 for normal galaxies (e.g. Stickel 
et al. 2000). This discrepancy becomes even higher, when 
considering the additional contribution of the atomic gas. 
Since it is unlikely that the CO-to H2 conversion factor for 
ULIRGs is so much different from that for normal galax- 
ies, a larger dust mass than M5^'°° is required for ULIRGs 
in order to match the canonical value. Such a suitably high 



Table 3. PAH 7.7 /im peak flux strength, 100 /xm and 850 jim 
continuum fluxes used for the assessment of flux ratios in 
Fig.^. Measured fluxes are in Jy. The MIR continuum around 
the PAH line was subtracted, also the extrapolated synchrotron 
continuum from the 850 /im flux. Uncertainties are less than 
30 %. For the ULIRG sample PAH fluxes are from Rigopoulou 
et al. (1999), 850 /im fluxes are from Rigopoulou et al. (1996), 
Lisenfeld et al. (2000) and from this work. For the sample of 
normal comparison galaxies PAH fluxes are from this work and 
850 fim fluxes are from Dunne et al. (2000). The 100 /im fluxes 
are aU from IRAS. 
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dust mass can be achieved with the combination of several 
dust components. 

In summary, we listed several items where the single black- 
body model appears to be insufficient or inconsistent. This can 
be solved by the application of multiple blackbodies. In order to 
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Fig. 5. Normalised 7.7 /im PAH flux versus opacity Tioo^m as 
formally derived from the single blackbody fits. Different sym- 
bols for each spectral type as in Fig.|. The arrows show the 
effect of extinction, if the PAH carriers are mixed with the FIR 
opaque dust. The filled circle and the thick vertical bar on the 
left-hand side show the range for normal galaxies "NGs" (see 
text). 



place several blackbodies underneath the Rayleigh- Jeans tail, 
/3 = 2 seems to provide a reasonable choice. 

5.1.2. Constraining the FIR opacity for the multiple blackbody 
model 

For a low FIR opacity and several cold to cool blackbodies 
with d = 2, the estimated total dust mass M,,, ,> is consid- 

' ' d(total) 

erably higher (on average by a factor of ten) than for the single 
optically thick blackbody (cf. Table This is because of the 
contribution from the additional cold component (M oc T^^ for 

(3 — 2 

P = 2). For the mm-subsample M^^^.^^^^^ lies in the range « 
10* - 10^ Mq. As mentioned in Sect. 4.5, the uncertainty for 



M 



13 = 2 
d(total) 



blackbodies is not unique. Nevertheless, when comparing the 
^d(7otai) ^'^^ '■^^ most reliable fits (marked by "H-" in Tab. ^ 
with the molecular gas mass derived from integrated CO lu- 
minosities (Tab.|^), the gas-to-dust mass ratio lies in the range 
between 30 and 300 with an average of 165 ± 120, close to the 
canonical value of about 150. The multiple blackbody model 
seems to yield a consistent total dust mass. 

Stricter constraints on the dust mass can be obtained, when 
comparing it with the extension of CO gas inferred from in- 
terferometric data (Downes & Solomon 1998). In order to 
keep the opacity, i.e. column density, low, the extent of the 
dust region must exceed a minimum size. As a check, we dis- 
tribute the dust mass associated with the bulk FIR emission. 



(3 — 2 

^d(FiR) (^^^^- r-^V evenly within a disk of minimal radius 
rr= 1, so that tiqq^„i = 1 (iV in Tab.|^. In each object the ra- 
dius r(Tioo^m = l) is smaller than the radius of the optical- 
NIR image sizes. With the possible exception of Arp220, 
the extent of low opacity dust is consistent with the size of 
the CO disk, where available (from Tables 3 & 4 in Downes 
& Solomon 1998): Mrk23 1 : 0'.'9, Mrk 273 : 0'.'9-3'.'l, 17208- 
0014 : 1'.'8, 23365H-3604 : 1'.'O. For noo;.™ = 0.3, corresponding 
to Ay ~ 50, the minimal radius r,- is about a factor two larger 
than listed in Tab. |[ and exceeds the CO disk sizes. If the dust 
and the molecular gas are well mixed, then the opacity of the 
bulk FIR dust component cannot be low (rioo/jm < 0.3), rather 
it must be "moderately" high (rJ^JJ^J^nj ~ 1, corresponding to 
Ay ~ 150) in many ULIRGs. Furthermore, with regard to our 
previous findings on the PAH/850 /im and PAH/100 /im ratios, 
which required relatively low extinction in the MIR (except for 
Arp220), we conclude that the bulk FIR dust component cannot 
contain the majority of the PAH carriers. The relative constancy 
of the PAH / 850 /im flux ratio with regard to normal galaxies is 
preserved, since the cold component dominates at 850 /im, as 
can be seen from the examples in Fig. |^. 

In summary, an emissivity exponent /3 = 2 works well for 
the ULIRGs, as in normal spirals and Markarian galaxies. In 
some ULIRGs the single component fits already give /3 = 2, and 
for those with /^fitted < 2 several blackbodies with /3 = 2 can be 
used. However, in order to match the dust mass and the size of 
the emitting regions with the constraints from the CO obser- 
vations, the opacity of the bulk FIR dust component might be 
moderately high (0.3 < t^^^^^ « 1). 

5.1.3. The proposed scenario: Cold cirrus dust and cool, partly 
opaque starburst dust 

The previous discussion suggests that the FIR-submm SEDs 
are composed of several (at least two) A^^ modified black- 
bodies with cool to cold temperatures. Actually, the dust might 
exhibit a continuous temperature range, but, when applying a 
formal model, within the uncertainties, a restriction to a few 
components works best. These dust components represent two 
stages: (1) the "cirrus" and (2) "starburst" stages. 

The nuclear regions are probably surrounded by clumpy 
"Super-Orion" complexes providing the heating power for the 
starburst dust. The gaps in between the cloud complexes and 



is quite large, because the decomposition into several 
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the outer regions with lower density and weaker interstellar ra- 
diation field host the cirrus-like cold dust. In this picture the 
PAH emission comes preferentially from the submm emitting 
dust. Some fraction of the FIR peak dust emission in the nu- 
cleus or in dense clouds can be opaque in the FIR. 

5.2. Starburst andAGN heated dust 

In this section we assess for each optical spectral class (see 
Tab. the typical shape of the IR-to-mm SEDs. While the 
SEDs of AGNs and SBs appear to be indistinguishable at FIR 
and submm wavelengths, they differ in the NIR-MIR. This sug- 
gests that the cool FIR emitting dust is not connected to the 
AGN, and that the AGN only powers the warm and hot dust. 
The difference at short IR wavelengths provides a NIR diagnos- 
tic tool based on the J-K colours, in order to reveal the presence 
of an AGN (with red colours, while SBs have blue colours). 
This tool is successful in the sense that if an AGN is seen by 
other diagnostics (optical or MIR spectroscopy or X-rays), then 
in most cases it is also seen via the NIR colours. 

5.2.1. Seyferts 

Our sample contains eight galaxies (with z < 0.3) which have 
optical Seyfert spectra, i.e. are sources known to house an 
AGN. Two are Seyfert Is and six are Seyfert 2s (see Tab. 

The SEDs of the Seyferts show a variety of shapes in the 
NIR and MIR range: 

• On the one hand, a power-law-like flux increase is seen in 
the Seyfert Is (Mrk231 and 15462-0450) and the Seyfert 2s 
05189-2524, 19254-7245 and Mrk463 (and probably also 
06206-6315). 

• On the other hand, Mrk273 (a LINER nucleus with an off- 
nucleus Seyfert 2 nebula, cf. Colina et al. (1999)) has a flat 
NIR flux plateau with a steep rise at about 10 /im. 

• 23389-6139, showing a damped power-law flux increase, 
appears to be in between these extremes. 

The power-law shape might be attributed to the central, rel- 
atively unobscured AGN. The red J-K colours might not be 
caused by extinction only, rather the hottest dust clouds at a 
temperature just below the evaporation temperature of the dust 
of about 1500K, correponding to a peak wavelength of about 
2 /im, might be seen. 

In a torus geometry the flat sources like 23389-6139 could 
be seen more edge-on, so that the hot AGN heated dust is more 
hidden. Then the flat NIR flux plateau - with blue colours - 
may be understood in terms of light scattered at bipolar cones. 
Alternatively, young hot stars in the host galaxies might con- 
tribute significantly to the NIR flux, shifting the colours to- 
wards the blue. Also, the AGN could be weak compared with 
the strength of the starbursts; perhaps Mrk273 belongs to this 
group. 

A somewhat indefinite case is the double source 06206- 
6315, which is classified by Due et al. (1997) as a Seyfert2, 
but as starburst via PAH diagnostics (Lutz, Veilleux & Genzel 



1999, Rigopoulou et al. 1999), although both spectra look quite 
noisy. The steep optical to NIR 2.2 /im flux increase argues in 
favor of a relatively unobscured AGN. 

In the FIR and submm range the Seyfert SEDs also ex- 
hibit some diversity, though only a relatively moderate one. Ex- 
amples are Mrk463 and 19254-7245, both with double nuclei 
(and similar power-law NIR flux increase): While Mrk463 only 
shows Httle FIR flux compared to the MIR flux, 19254-7245 
has a strong FIR bump and even a considerable amount of mm 
flux due to cold dust. A similar difference is found for Mrk23 1 
and 05189-2524, both with a "single" nucleus, i.e. probably 
mergers in an advanced state: Both show a strong FIR bump, 
but Mrk23 1 has a steeper Rayleigh- Jeans tail and less cold dust 
than 05 1 89-2524. Thus, the host galaxies of AGNs may or may 
not have cold dust. 

5.2.2. LINERs and HII/SBs 

Our sample contains 11 LINERs (Low lonisation Nuclear 
Emission Regions) and 9 HII/SBs. In general, the SEDs look 
very similar for both optical spectral types. Therefore, we dis- 
cuss them together here. 

The SEDs appear quite homogeneous from the NIR to the 
FIR, with some diversity in the submm range: All SEDs show 
a relatively flat NIR flux plateau with a faint 2 /im bump, PAH 
emission and/or Si 9.7 /im absorption features, then a steep 
rise in flux at about 10 /im, i.e. a 10 /im knee. The flux peaks 
in the FIR at about 60-1 00 /im. The only variety appears in 
the submm range: Some sources show a steep Rayleigh-Jeans 
tail, others exhibit a flatter tail providing evidence for addi- 
tional cold dust. Typical sources with and without cold dust are 
23230-6926 and 23365H-3604, respectively, among the LIN- 
ERs, and 17208-0014 and ES0286-19, respectively, among the 
HII/SBs. 

The homogeneity and extreme similarity of the NIR to FIR 
SEDs suggests that the dust in LINERs and HII/SBs has sim- 
ilar properties, in particular concerning the spatial distribution 
and heating mechanism. In fact, this similarity places the LIN- 
ERs closer to HII/SB than AGN dominated ULIRGs. A similar 
conclusion (that LINERs are not AGN dominated) was drawn 
from the PAH diagnostics (Lutz, Veilleux & Genzel 1999). 

5.2.3. The SED shapes and J-K colours as diagnostic tool 

The NIR-MIR power-law flux increase for the AGN-type 
ULIRGs and the flat NIR flux plateau with a 10 /im knee for the 
SB-type ULIRGs suggests the utilization of these two different 
SED shapes for a classification scheme. Earlier investigations 
(e.g. Hill et al. 1988) showed that red NIR colours provide a 
high probability of finding an AGN among moderately lumi- 
nous IRAS galaxies. As a quantitative characterisation mea- 
sure for the SED shapes we use the J-K colours: J - K > 2 mag 
identifies power-law shape SEDs; this value corresponds to a 
2.2/im/ 1.25/im flux ratio of about 2.34. 

Although longer MIR wavelengths are less sensitive to ex- 
tinction, this range suffers from possible confusion by strong 
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PAH emission and 9.7 /im silicate absorption which compli- 
cates the analysis; detailed work has been done by several au- 
thors, e.g. Laurent et al. (2000), Imanishi & Dudley (2000) and 
Tran et al. (2001). 

As shown in Fig.^ the J-K colours (2.2/im/ 1.25/im flux 
ratios) provide a much clearer separation of AGN- and SB- 
ULIRGs than the MIR colours (25^m/60^m) used to identify 
warm AGN-type objects (de Grijp et al. 1985), to say nothing 
of the FIR colours (120/J,m/ 200/im). 

In Fig. ^ we show in histograms the correspondence of this 
NIR photometry classification scheme with the spectroscopic 
optical emission line and the MIR PAH diagnostic classifica- 
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shape classification (via the J-K colour) with the optical emis- 
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tion as compiled in Table Q With respect to the optical classi- 
fication, we find that the NIR SED shape confirms for 22 out 
of 22 (100%) the optical classification as a HII/LINER object 
and for 6 out of 8 (75 %) the classification as Seyfert 1/2. With 
respect to the MIR PAH diagnostics, we find that the NIR SED 
shape confirms for 12 out of 14 (86%) the PAH classification 
as SB and for 4 out of 5 (80 %) the PAH classification as AGN. 

These classifications are further confirmed by comparing 
them with the hard X-ray classifications (Table^. Thus, the J- 
K colours represent an attractive tool to explore the nature of 
ULIRGs and other dust rich IR galaxies. This is of particular 
advantage, if the galaxies are too distant or too faint for spec- 
troscopy. Then the intrinsic J-K colours can provide constraints 
on the nature of cosmologically interesting ULIRGs. 

All these methods have some limitations, of course, and 
in some cases provide indefinite classifications. For example, 
the J-K colour does not reveal any AGN in Mrk273 which is 
optically classified as Seyfert 2 and has a PAH/continuum flux 
ratio of 1.9 favouring the starburst dominance (Rigopoulou et 
al. 1999). Also X-ray observations with Beppo-Sax unveiled 
an AGN in NGC6240 (Vignati et al. 1999), which is optically 
classified as a LINER, SB dominated according to PAH diag- 
nostics, and shows no AGN-typical power-law NIR-MIR flux 
increase (except for the 1" nuclear region). 



U. Klaas et al.: Ultra-luminous IR galaxies 



19 



5.2.4. The nature of the 10 /im knee 

Pure starburst objects have a flat NIR flux plateau followed by 
a steep rise longwards of 10 /im. Part of this pronounced 10 /im 
knee could be due to the 9.7 /im silicate absorption. In contrast, 
galactic compact HII regions housing 05-09 stars show SEDs 
with a steep power-law-like rise. As for AGNs this rise already 
starts at about 1 /im (Chini et al. 1987). Also, the prototype 
starburst galaxy M82 has such a SED shape (e.g. Thuma et al. 
2000). Thus, a power-law SED could, in principle, also occur in 
ULIRGs with starburst classification. But this is not observed 
for our sample. 

HST NICMOS images (Scoville et al. 2000) reveal that 
the central 1" region is redder than the outer regions. 
Fig.|l] shows that for those 9 sources in common with 
our sample (05189-2524, 12112H-0305, Mrk273, 14348-1447, 
15250H-3609, Arp220, NGC6240, 17208-0014 and 22491- 
1808) most of the nuclear SEDs (marked with "x" symbols) 
exhibit - even for LINERs and HIIs - a power-law-like flux 
increase in the NIR. J-K colours of these nuclear regions are 
similar to the colours in larger apertures found for the AGN- 
type ULIRGs. For 05189-2524 the NICMOS fluxes (hardly to 
recognize in Fig. |l]) are practically identical with the ground 
based fluxes showing that this source is completely dominated 
by the central point-like source. For the others the much higher 
and bluer ground-based fluxes suggest that the sources are ex- 
tended on the 5 - 10" scale, in particular for 17208-0014 and 
22491-1808 for which already the NICMOS images indicate 
that the SEDs are composed of many different components. 

With regard to the Balmer decrement LINERs and HII- 
starburst galaxies have low Eb-v (Tab.^, as found for 
Seyfert Is. The highest Eb-v are those for Seyfert2s, indicat- 
ing that their emission regions are the most deeply embedded 
ones. A deeply embedded HII region should suffer from extinc- 
tion in the same way. However, deeply embedded young stars 
will probably be outshone by those located closer to the outer 
surface of the galaxy and will not dominate the NIR colours. 
In conclusion, the relatively flat blue NIR flux plateau comes 
from the practically unobscured outer regions. The strong flux 
rise longwards of 10 /im is caused by more deeply embedded 
luminous star forming regions. This naturally explains the ap- 
pearance of the 10 /im knee. 



5.2.5. Does the AGN power the FIR emission of ULIRGs? 

It is now widely accepted that a powerful (and not obscured) 
AGN creates a strong NIR-MIR emission (e.g. Sanders et al. 
1988a, 1988b, Pier & KroUk 1992, 1993, Rowan-Robinson 
1995, Haas et al. 2000a). However, it is still a matter of de- 
bate whether the FIR/submm emission in ULIRGs (as well 
as quasars) is also mainly powered by the AGN (e.g. Sanders 
1999) or by circumnuclear SBs (e.g. Rowan-Robinson 1995, 
Genzeletal. 1998). 

The power-law flux increase of Mrk23 1 appears similar to 
that of the quasars PG0050H-124 and PG 1613H-658 (Haas et 
al. 2000a). For these quasars it is not clear whether and in 



which wavelength range a starburst contributes to their power- 
law SEDs: Energetically, the AGN radiation would be suffi- 
cient to heat the dust, and since pure starbursts have different 
SED shapes (with knees, as shown above), the superposition of 
a significant starburst component would probably destroy the 
smooth power-law shape. For Mrk231, however, the FIR lu- 
minosity exceeds the MIR- and the (dereddened) UV-optical- 
luminosity (Downes & Solomon 1998). Downes & Solomon 
(1998) found two rotating gas disks with radii of 460 and 1 150 
pc, respectively, and concluded that the outer one is starburst 
dominated, since it is too extended and FIR-luminous to in- 
tercept sufficient power from the central region. Although the 
AGN is clearly visible, it dominates neither the FIR luminos- 
ity nor the total energy output. However, the PAH/continuum 
diagnostics indicates a dominating AGN in Mrk23 1 (Genzel et 
al. 1998, Lutz, Veilleux & Genzel 1999). The reason for this 
might be that the unobscured AGN continuum diminishes the 
relative PAH strength, mimicking a weak starburst and overem- 
phasising the role of the AGN. Thus, it is necessary to consider 
not only the ratio PAH/continuum, but also the absolute PAH 
and continuum levels, both of which may be affected by differ- 
ent amounts of extinction due to different sites of origin in the 
galaxies. For the diagnostics of moderately luminous Seyfert Is 
and Seyfert 2s this has akeady been stressed by Clavel et al. 
(2000). 

As shown in Fig.||, the ULIRGs with and without AGN 
signatures cannot be distinguished via their FIR colours 
(Fi20/^m/F2oo^m)- This suggests that either the FIR/submm 
emission is largely independent of the presence of an AGN, 
or that every ULIRG not classified as AGN contains a hid- 
den AGN. The relative high number of ULIRGs with signs 
of strong starbursts versus ULIRGs with AGN signatures in 
our sample (22:10 from optical emission line diagnostics, 13:5 
from PAH diagnostics) would argue in favour of the first al- 
ternative. Since in Sect. 5.1.2 we concluded that the bulk FIR 
component might be moderately opaque (Tj^o"o%n ~ 1 ) even IR 
diagnostics, in particular the PAH / 850/im flux ratio, have lim- 
ited meaningfullness in this context. An exception seems to be 
Arp 220 where the dust geometry leads to a strong deficiency 
of PAH strength versus 850 /im flux and dereddening would 
yield a quasar-like IR continuum (Haas et al. 2001). Again, 
this is only 1 out of 22 sources of the submm ULIRG subsam- 
ple showing this behaviour Therefore, we see some evidence 
from our sample that the role of the AGN for powering the 
FIR/submm emission in nearby ULIRGs is negligible. 

5.3. Cold dust and search for evolutionary trends 

ULIRGs may well have considerable amounts of cold dust at 
10-30 K, a temperature typical for dust in less active spiral 
galaxies. If the cold dust (from the parent galaxies) is continu- 
ously heated during the merging process, then we would expect 
more advanced mergers to show a relatively larger amount of 
warm dust than mergers in a beginning phase. To check this 
hypothesis, we use the morphological appearance and the sep- 
aration of the galaxy nuclei as a measure for the merger state. 
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This is a simplification, since the encounter might not lead to 
a monotonous approach of the two galaxies. According to sim- 
ulations by Dubinski, Mihos & Hernquist (1999) the merging 
process contains repeated approaches with semi-elastical col- 
lisions followed by a drift apart. Nevertheless, we divided the 
sources into (1) single and (2) double (see Tab. We use the 
MIR, FIR and submm luminosities as a measure for the amount 
of warm, cool and cold dust, respectively. Fig. || does not show 
any trend of M1R-, FIR- and submm-luminosity ratios with 
the simple morphological classifications. Within our sample, 
which covers only a small range in bolometric luminosity, the 
relative amounts of warm, cool and cold dust are not correlated 
with the merger state inferred from the morphology. 

A similar independence was found between the strength 
of the PAH-feature/MIR-continuum and the (projected) sepa- 
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ration of the galaxy nuclei (Lutz et al. 1998, Rigopoulou et 
al. 1999). Either the (projected) separation of the nuclei is not 
well suited to measure the progress of the merging process, or 
the dust heating does not steadily increase during the merging. 
It seems likely that while the merging nuclei approach each 
other relatively slowly and/or repeatedly, the dust clouds un- 
dergo several phases of compression, turning into cold proto- 
starburst clouds which subsequently undergo heating by star- 
bursts (and by an AGN, if present). These distinct phases could 
occur simultaneously in a merger, but at different locations. 
Photometry of the entire galaxy then shows a mixture of dust 
complexes at different phases. This is actually revealed by spa- 
tially resolved submm observations of the colliding IR lumi- 
nous (not yet ultra-luminous) galaxy pair NGC 4038/39 (Haas 
et al. 2000b). 

Irrespective of the morphology, the F25/F6Q versus 
Fi5o/F85o two-colour diagram (Fig.^ illustrates the distribu- 
tion of our ULIRG mm-subsample. Schematically, the cold 
cirrus-like dust is represented by the 850 /im emission, the 
cool SB dust by the 150 and 60 fim emission, and the warm 
AGN dust by the 25 /im emission. In this diagram trends 
are recognizable: strong AGNs are located in the upper half 
(F25/F60 JlO.2), SBs in the lower half. The galaxies with and 
without cold cirrus-like dust lie towards the left and right side, 
respectively. This diagram illustrates that despite similar lumi- 
nosities the nearby ULIRG sample exhibits quite some diver- 
sity. Whether or not this is due to evolution is still a puzzle. 

Finally, it is interesting to note that the four ULIRGs at 
medium redshift (z w 0.3) have similar dust temperatures as 
the low redshift ULIRGs, but exhibit the highest FIR lumi- 
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nosities among the sample. While their MIR SEDs and their 
spectral types can be determined soon, the question of whether 
they also contain cold dust will have to be left to future obser- 
vations with the Herschel Space Observatory or the Atacama 
Large Millimetre Array. 

6. Summary 

Infrared to millimetre spectral energy distributions (SEDs) 
have been obtained for 41 bright ultra-luminous infrared galax- 
ies (ULIRGs). These are the most up-to-date and detailed pho- 
tometric templates of the nearby ULIRGs, which are the fun- 
damental ingredient for cosmological studies. Arp220, which 
has been considered so far as a the archetypal ULIRG, is one 
of the most peculiar objects in our sample. 

The SED maxima lie between 60 and 100 /im, and for those 
22 sources with submm detections or upper limits the slopes of 
the Rayleigh-Jeans tails can be well constrained. The FIR and 
submm parts of the SEDs between 60 and 1 300 ^m can be fitted 
in two ways: 

Firstly, with a single modified blackbody yielding large 
ranges for both the opacity (0.1 < rioo/^m < 5) and the emissiv- 
ity (1.2 < /3 < 2.2). The resulting temperatures range from 50 
to 70 K. However, there are several reasons against the physical 
relevance of one single dust component: 

1) The optical-NIR morphology is generally disturbed, sug- 
gesting the presence of diffuse outer dust layers in addition 
to a dense central dust body. 

2) It is difficult to understand why (3 is significantly < 2 in 
some ULIRGs, but /3 « 2 for spiral and Markarian galaxies. 

3) The normalised strength of the 7.7 /im PAH features shows 
no correlation with the extinction derived formally from 

4) The derived dust masses lead to unusually high gas-to-dust 
mass ratios from 500 to 2300. 

Secondly, the FIR-submm SEDs can be decomposed into 
two or more A^^ modified blackbody spectra representing cool 
(30-50 K) and cold (10-30K) dust components. While the cool 
component should be moderately opaque (0.3 ^ t}}5oJ^j^ wl)to 
match size constraints inferred from interferometric CO obser- 
vations, the cold component should be transparent in order to 
conform to the observed PAH/100 /im and PAH/850 /im ratios. 
Remarkably, this suggests that the bulk FIR dust component 
cannot contain the majority of the PAH carriers, rather they are 
more closely coupled with the submm emitting dust. 

This second way to fit the SEDs provides evidence favour- 
ing two dust stages: 

1) the cold dust in more quiet regions, which is mainly respon- 
sible for the submm emission, and 

2) the cool dust powered by starbursts (SB) which dominates 
the FIR and may be partly optically thick even in the FIR. 

The dust masses derived from such a multiple component 
model yield gas-to-dust mass ratios in the range between 30 



and 300, hence close to the galactic value of 150. ULIRGs clas- 
sified as either AGNs or SBs appear indistinguishable at FIR 
and submm wavelengths. We suppose that the cool FIR emit- 
ting dust is not related to the AGN, in particular since SB-type 
ULIRGs seem to be more than twice as frequent in our sample 
as AGN-type ULIRGs. 

In the NIR-MIR the SEDs reveal two basic shapes: 

1) A flat NIR-MIR flux plateau followed by a steep rise at 
about 10 /im. All objects optically classified as LINERs 
and HII/SBs exhibit this SED shape, but so do some of the 
Seyferts. 

2) A power-law like flux increase starting in the NIR and 
peaking in the FIR or even MIR. This SED shape is only 
seen in Seyferts. 

The NIR-MIR power-law-like emission is due to warm and hot 
dust powered by the AGN, and this is a third dust stage re- 
stricted to AGN-type ULIRGs. From this finding, a NIR pho- 
tometry diagnostic tool is established to reveal the presence of 
an AGN on the basis of the J-K colours (which has red colours, 
while SBs have blue colours). This tool is successful in the 
sense that if an AGN is seen by other diagnostics (optical or 
MIR spectroscopy), in most cases (6 out of 8) it is also seen via 
the NIR colours. 

The detection or non-detection of cold dust could reflect 
different evolutionary states of the ULIRGs associated with 
the merger state. However, the current means to determine the 
merger state, like the projected distance of the nuclei, appears 
to be insufficient to verify evolutionary trends. 

Appendix A: Near-Infrared images 

The J- and K-band images look very similar, therefore only the 



K-band images are shown (Fig. Most of the emission is 



concentrated, but also irregular/extended tails are present. The 
pointlike blue source 10" southwest of 00262H-4251, and the 
two ones located about 6" north and south of 18090H-0130 are 
probably stars (and are excluded from the photometry). The 
photometry derived from the images is listed in Tab.[l|. 

Appendix B: Optical spectra 

Fig. B.l shows optical spectra for two sources, obtained with 
the TWIN spectrograph at the 3.5 m telescope on Calar Alto, 
Spain. Both sources have a LINER spectrum according to the 
line diagnostics by Veilleux & Osterbrock (1987). Now all 
sources of our mm-subsample are spectroscopically classified 
(see Tab. |). 
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Fig. B.l. Spectra of 20046-0623 (top) and 21504-0628 (bottom), blue (left) and red (right) channel of the TWIN spectrograph. 
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Table 4. Morphologies and spectroscopic classifications. The 
morphologies are from optical and NIR images as referenced. 
The extinction Eb-v has been compiled from H^/H/j ra- 
tios. The optical spectral types distinguish between Seyfert 
1, Seyfert 2, LINER (galaxies with Low lonisation Nuclear 
Emission Regions), and Hll/starburst types, as classified in 
the respective references. The MIR spectral types distinguish 
between AGN = Seyfert 1 h- Seyfert 2 and SB = LINER h-HII 
types (Lutz, Veilleux & Genzel 1999). The NIR classifications 
are derived from the J-K colours (AGN, SB, SB/n = SB for to- 
tal galaxy and red 1" nucleus, SB/f = SB for total galaxy and 1" 
nucleus). The X-ray evidence for AGN is taken from Risahti et 
al. (2000). 



Name Morph Eb v Optical MIR NIR X-ray 





[mag] 


Sp-type 


type 


'^^^ 


AGN 


00199-7426 double 




LINER? 




SB 




00262-^4251 merger '"^•* 








SB* 




00406-3127''i merger """^ 




Seyf 2 








03068-5346 double?"^'' 












03158+4227 compact"^ 










No 


03538-6432''^ compact"^'^ 












04232-fl436 merger * 








AGN* 




05189-2524 merger™^ 


2.03 


Seyf 2 


AGN° 


AGN 


Yes 


06035-7102 double 


1.34 


HII/SB =1 


SB 


SB 




06206-6315 double 


2.06 


Seyf 2 "1 


SB 


AGN 




12112+0305 double 


0.92 


LINER °^ 


SB 


SB/n 




Mrk 231 merger 


0.70 


Seyf 1 


AGN 


AGN 


Yes 


Mrk 273 merger 


1.22 


Sy2 "'iLm'^ 


AGN 


SB/n 


No 


Mrk 463 double"*" 


0.60 


Seyf 2^'' 


AGN 


AGN 


Yes 


14348-1447 meiger ""^ 


1.15 


LINER 


SB 


SB/n 




14378-3651 meiger ""^ 


1.26 


LINER ^1 




SB 




15245+1019 merger ""^ 












15250+3609 ring gal. 


0.76 


LINER 


SB 


SB 




Arp 220 merger 


1.05 


LINER 


SB 


SB/n 


No 


15462-0450 merger 


0.60 


Seyfl 




AGN 




16090-0139 merger "•^■^ 


1.55 


LINER "^'^ 




SB 




NGC 6240 double 


1.65 


LINER 


SB 


SB/n 


Yes 


17208-0014 merger "^'^ 


1.75 


HII/SB 


SB 


SB/n 


No 


17463+5806''i compact""^'^ 












18090+0130 merger * 








SB* 




18470+3234 double"^'* 


1.13 


HIl/SB =^ 




SB* 




19254-7245 double "^-^ 


2.92 


Seyf 2 =1 


AGN 


AGN 


Yes 


19458+0944 double 












20046-0623 double "^■^•* 




LINER * 




SB 




20087-0308 merger "^'^ 


1.43 


LINER 




SB 




20100-4156 double 


1.12 


HII/SB 


SB 


SB 




20414-1651 merger? "i'^ 


1.31 


LINER '1 








ESO 286-19 merger 


1.07 


LINER ^1 


SB 


SB 


No 


21130-4446 merger"*^ 


0.50 


Hn/SB '1 




SB 




21504-0628 meiger"*^'* 




LINER* 




SB* 




22491-1808 meiger ""^ 


0.81 


Hn/SB 


SB 


SB/f 




ESO 148-2 merger ""^ 


1.01 


Hn/SB 


SB 


SB 


No 


23230-6926 merger ""^ 


1.49 


LINER 




SB 




23365+3604 merger'"'''* 


0.88 


LINER*'' 


SB 


SB* 




23389-6139 double "^^ 


3.06 


Seyf 2*1 




SB 




23515-3127''i compact"®^ 




Seyf 2*'' 








'•^ z>0.3 

* this work, see Appendix 
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